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In a recent letter Q (referred to as I below), Machida 
et al. made the exciting claim that in a one- dimensional 
(ID) trapped gas of fcrmions with repulsive interactions a 
supcrfluid phase appears around the Mott-insulator (MI) 
at the center of the trap (COT). Their claim is based on 
a negative binding energy (Eb), and a large weight for 
a singlet formed by particles located at opposite sides of 
the MI. We show here that the observed effects are not 
related to superfluidity. 

After a MI forms at large U, two particles with oppo- 
site spins added to the trap prefer to sit beyond the two 
ends of the MI phase in order to avoid double occupancy. 
Hence, the large weight of the singlet [Eq. (3)] in Fig. 4 
of I can be understood to be a simple consequence of the 
density distribution and the antiferromagnetic character 
of the MI state, i.e., it does not signal superfluidity. 

We then focus on the origin of the negative Eb ob- 
served in I. Most of the results in I exhibit a non-zero 
density at the borders of the trap, i.e, they depend on 
the boundary conditions. We thus recalculated two cases 
depicted in Fig. 1(a) of I, keeping the same curvature of 
the trap and increasing the system size to TV = 20. This 
ensures zero density at the borders, as it should be for 
confined systems. We used quantum Monte Carlo (QMC) 
simulations 0, , density-matrix renormalization group 
(DMRG) 0, and exact diagonalization (ED). 

Figure ^ shows that a negative Eb appears for large 
values of U/t, at the point where the MI sets in the COT. 
Beyond this point adding more particles to the system 
causes the MI to increase in size. This is in contrast to 
the systems without a trap, where adding more particles 
to the half-filled case causes the MI to disappear. Both 
negative and positive Eb arise in the latter doped case 
for different boundary conditions 

The dashed line in Fig. 1 corresponds to ED results of 
the ID Hubbard model without trap, at half-filling, and 
open boundary conditions (OBC). Here Eb is calculated 
by adding a site when adding a particle, for OBC, in or- 
der to simulate the MI in the middle of the trap without 
the metallic wings. The results obtained are practically 
indistinguishable from the Eb obtained for trapped sys- 
tems after the MI appears in the COT. Therefore, the 
negative Eb is due to the MI region, and does not sig- 
nal superfluidity in the wings of the MI. Moreover, in 
the inset of Fig. ^ we show that the negative Eb in the 



FIG. 1: (color online). QMC (points) and DMRG (continuous 
lines) results for E b (W) = E(5, 5) + E(6, 6) - 2E(5, 6) vs U/t 
for: V = 7.5t (V,6)> and V = 9.5t (OA) (see I), when the 
center of the trap is located in the middle of two lattice points 
(\7,0), and in a lattice point (Q,A). The dashed line shows 
ED results for Et(10) in open MI systems (see text). The 
inset shows ED results for Eq — Eg + ANf vs TV/, in open 
MI systems with U/t = 8. Here, we choose a nonzero value 
of A = 0.327 to stress the even(x)-odd(+) effect, without 
affecting the actual value of Eb < it causes. 


MI is due to an even-odd effect. There we have plotted 
the ground-state energy Eq for MI systems without the 
trap, and OBC, vs the number of particles (TV/). The 
even-odd effect is evident, and becomes smaller with in- 
creasing system size. Additionally, consistent with the 
results above, we find that: (i) Displacing the COT from 
the middle of two lattice points, as selected in I, leads to 
positive values of Eb- Results for the COT on a lattice 
point are also shown in Fig. ^ (ii) In the trap, similarly 
to the OBC case, the negative Eb — > almost linearly 
with increasing system size, i.e, it is a finite size effect. 
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